Regularization of Diagrammatic Series with Zero Convergence Radius 
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The divergence of perturbative expansions for vast majority of macroscopic systems, which follows 
from Dyson's collapse argument, prevents Feynman's diagrammatic technique from being directly 
used for controllable studies of strongly interacting systems. We show how the problem of di- 
vergence can be solved by replacing the original model with a convergent sequence of successive 
approximations which have a convergent perturbative series. As a prototypical model, we consider 
the zero-dimensional \Tp\'^ theory. 
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Almost sixty years ago, Dyson provided a physical ar- 
gument why power-series expansions in quantum elec- 
trodynamics are divergent even after proper charge and 
mass renormalization [1]. The argument is straightfor- 
wardly generalizable to a wide range of classical and 
quantum statistical problems and has far-reaching con- 
sequences for any perturbative approach J2j. It goes as 
follows: The interaction potential between two electrons 
with charge e is of the form e^/ri2, with ri2 the distance 
between them. When performing a series expansion in 
{i.e., the interaction strength) a physical quantity of 
interest can be expressed as 

F{e^) = ao+a2e'^ + aie* + ... (1) 

If the series have finite convergence radius ^ then F{e^) 
is an analytic function at e — 0, implying that for 
sufficiently small values of e, F(— e^) is a well-behaved 
analytic function. However, F{—e^) corresponds to a 
fictitious world with purely imaginary particle charges 
where the interaction potential is of the form — e^/ri2, 
and is thus attractive. In this fictitious world, the 
vacuum state would be unstable against production of 
an infinite number of electron-positron pairs, each put in 
a separate region of space since the gain in the negative 
Coulomb energy is larger than the increase in the kinetic 
energy. Hence, Dyson concluded that F{e^) cannot be 
analytic around e = and that the convergence radius 
must be zero, ^ = 0. 

The 1-01^ theory is associated with the following parti- 
tion function in d dimensions: 




describing the statistics of the classical complex-valued 
field ip{r). This effective field theory is often used to de- 
scribe critical behavior of the superfluid phase transitions 
in interacting Bose systems. The model ([2| can be solved 
numerically very efficiently by a number of methods, but 
none of them is based on an expansion in A. A successful 



regularization of the perturbative expansion in A for the 
IV' I'' theory would change the status of the standard dia- 
grammatic technique to that of a systematic method for 
accurate studies of strongly correlated systems, and pro- 
vide the basis for numerical treatments in the framework 
of (bold) Diagrammatic Monte Carlo [3] . Since Taylor ex- 
pansion in A is subject to the Dyson's collapse argument, 
it is clear that some regularization of the theory is re- 
quired by which we mean any sequence of field-theoretical 
approximations to the original model, controlled by reg- 
ulator N, leading to a diagrammatic series with (at least) 
a finite convergence radius. The correct physical result is 
obtained after extrapolation of N to infinity. As an ex- 
ample of such a regularization, we mention a continuous- 
space (finite-temperature) system of fermions with trun- 
cated single-particle momenta. Truncation sets the limit 
on the maximal density of the system, thereby saving the 
problem from Dyson's collapse and allowing one to em- 
ploy the diagrammatic treatment up to high-order (see, 
£■9-, [l])- While we do not exclude the possibility that 
in specific cases the regularization procedure might ulti- 
mately be reduced to a re-summation scheme, we real- 
ize that a generic regularization protocol, in view of the 
mathematical ambiguity of restoring a function from an 
asymptotic series with zero convergence radius, should 
derive from the series provenance, i.e. explicitly depend 
on the form of the original theory. We stress that we 
seek a regularization which preserves the structure of di- 
agrammatic expansions. 

In this Letter, we introduce a class of regularization 
techniques that introduce counter-terms to the diagram- 
matic series. In this approach, the so-called "sign bless- 
ing" , i. e. sign alternation of the same-order diagrams 
(cf. the fermionic case [5]) of a regularized series — as 
contrasted to the sign-definiteness of the same-order di- 
agrams of the original theory — is responsible for mutual 
cancelation of contributions from the factorial number of 
diagrams of a given order and series convergence. 

Focusing solely on the convergence properties, it is 
instructive [2J to reduce the theory ^ to its zero- 
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dimensional extreme. In this case, the field ipi^) is re- 
placed by a complex number ip and the partition function 
reduces to the integral {x = IV'P) 




In order to understand the convergence properties of the 
perturbative approximations in A for the 1-0 1** theory, it 
suffices to do that for the integral /(A). 

We start with identifying the problem for the standard 
Taylor series representation 

e-^^ = hm y x''' , (4) 

Ar-)-oo^ fc! 
fc=0 

which is that, reversing the order of taking the limit 
N oo and integrating over x, one obtains a series di- 
verging as~ A'"'(2fc)!/A:!. This is a hallmark for an asymp- 
totic series with zero convergence radius. An equivalent 
observation is that for any finite N the integral is domi- 
nated by values of x where the finite sum is not providing 
an accurate description of the original exponent ^ . Us- 
ing a different representation for the exponential function 

e-^"' = lim Zn{x)= lim (1-^) , (5) 

also fails for exactly the same reason. It formally pro- 
duces a sequence of integrals InW, by replacing e~^^ — >■ 
Z]^{x), which takes us away from the original problem 
because the dominant contribution to /at (A) comes from 
X ^ N, as contrasted to the scale a; < 1 of the original 
integral. Equation Q can be transformed into Eq. ^ by 
applying a resummation technique to the original Taylor 
series, 

(l-^)'^-E^^--^V(fc,iV), (6) 
^ ^ fc=0 

where the resummation function f{k,N) — k\C^/N^ is 
such that f{k <C N) « 1 and decreasing fast for k N. 
Formally, / is supposed to suppress the leading diver- 
gence of the original series which is not the case for 
Eq. ^. One may seek other functions /(fc, N) rendering 
the extrapolation procedure N ^ oo meaningful but as 
far as we know the solution was not found yet, nor is it 
clear that it exists. However, as we prove below, the solu- 
tion does exist when the series are generalized to include 
additional counter-terms; namely, within the 

A/— ^oo ^ — ^ 

fc=0 

representation, where /(fc, N) is restricting summation 
over fc to some polynomial of finite order which scales 
with N. We require that the expansion is in terms of 



integer powers of x since otherwise the construction of 
Feynman's diagrams (which is our ultimate goal) be- 
comes problematic. To generate the standard Feynman 
diagrammatic expansion based on Wick's theorem one 
has to introduce the exponential 

Zjv(a;) = e'"[^"(")] , (8) 

and expand the logarithm in powers of x. 

Our solution is based on designing an alternative rep- 
resentation of the exponential function wich has to sat- 
isfy several requirements when expanded in powers of 
the coupling parameter. A successful perturbative ap- 
proach has to satisfy the conditions of (i) convergence of 
the A-expansion of /at (A) for all A > at any fixed TV, 
and (ii) meaningful extrapolation to the N ^ oo limit 
to ensure that limjv-i-oo -^w (A) — /(A). Given that the 
theory ^ features non-analyticity at the second-order 
phase transition point, one might think that the condi- 
tion (i) should be weakened to sufficiently small A's. Nev- 
ertheless, we consciously require convergence everywhere 
to guarantee that the divergence of the A-expansion of 
the field-theoretical counterpart of /at (A) is exclusively 
due to the phase transition in the macroscopic limit; the 
convergence radius remaining infinite for any finite-size 
system (qualitatively analogous to our zero-dimensional 
model). 

Let us introduce the parameterization 

= u=(A/iV)i/", (9) 

with TO an integer number, and note that the most se- 
vere problem with Eq. ([s]) can be overcome by seeking a 
solution in the form 

e-^-^ = hm [fiy)f , (10) 

A*— >oo 

with the function /(.) being bounded, < /(.) < 1, along 
the real axis of y, and having an expansion 

/(y) = 1 - + . . . . (11) 

The terms that are of a higher exponent than m serve as 
counter terms compared to Eq. ([5|. These requirements 
guarantee that the series of integrals InW converge to 
the right answer, but this does not solve yet the problem 
with perturbative expansions unless we also require that 
the integral 

poo 

/jv(A)= / dxe-^[f{y{x))f (12) 
Jo 

is convergent for any complex value of u [see Eq. Q], 
which is the necessary condition for /at (A) to be expand- 
able in the convergent Taylor series in powers of u. It is 
this condition which forces one to consider to > 1. 

It turns out that there are infinitely many functions 
corresponding to to = 4 (to = 4 is in practice the obvious 
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choice), all satisfying the above- formulated requirements. 
Here we present one of them which is based on Bessel 
functions J„: 



f{y) = .M2z) + 2J2{2z) + -M2z) , 



(13) 



where z — (72)^/"* y and the explicit values of coefficients 
follow from the requirement (11). With regards to all 
conditions we find; 

- For any positive u the condition </(.)< 1 means 
that all integrals ( 12 1 are dominated by the a; < 1 region 



and thus guaranteeing that the condition (ii) is satisfied. 
- [f{y)]^ is an entire function with infinite convergence 
radius. 



- The integrals ( 12 ) converge for any complex u since the 
strongest divergence at x — ?> cxd (occurring at imaginary 
u) is exponential with the exponent oc and thus not 
dangerous in view of the e~^ factor; the condition (i) is 
thus satisfied. 

- the Taylor expansion of f{y) is in even powers of the 
argument. Hence, the requirement of having only integer 
powers of x is also satisfied since y^ = u^x. 
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FIG. 2: Plotted are integrals Ip^M = J dxe ^ Zp^m{x) where 
Zp^N is the Taylor expansion of [f{y)\^ up to order , for 
A = 0.5. In this example N — 2 which guarantees that the 
converged answer is accurate down to one percent accuracy. 
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FIG. 1: (color online) Convergence properties of the sequence 
/iv(A) with f{y) defined in Eq. ([l3|. The result for /(A) is 
shown with circles. The convergence with A*' for A = 5 is 
shown in the inset. 



FIG. 3: Integrals Jp.jv = / da;e~^Zp.jv(a;) based on the same 
series as in Fig.[2]resummed using the near-Gaussian function, 
see text. The elimination of the severe sign-problem and the 
radical improvement of convergence properties suggests that 
polynomial size and optimization should be explored. 



The convergence with N of the sequence of integrals 
( 12 1 is illustrated in Fig. [I] Remarkably, iV = 4 already 
produces results which are accurate at the 1% level up 
to A = 5. The ultimate convergence to the exact solu- 
tion is a polynomial function of N^^^^ and can be easily 
extrapolated from finite N. A peculiar feature of the 
arising theory is that perturbative expansion is in terms 
of u which is a fractional power of A. 

In order to comply with Feynman's diagrammatic rules 
we need to study integrals over x for Taylor expansions 



Zn{x) 



Tl'kLo^^^k done up to order p. The ir 



plications for that are shown in Fig. [2] where we plot 
the final result when Z]\[{x) is replaced with Zp^^^^x) — 
J2k=o ^'^^k- At first, the impression is that a severe sign- 
problem and a very high-order expansion is the price for 
obtaining a set of approximations using series with in- 
finite convergence radius. When we take stronger A or 
larger TV, the fiuctuations that need to cancel against 
each other become larger, and the minimum Taylor ex- 
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pansion order increases. However, when the original se- 
ries are further modified using resummation technique 
based on the near-Gaussian function 

p 

Zp,Ar(a;)^Zp,jv(a;)=^a;*'-Cfce-'^°/P, s = 2.05, (14) 

fc=0 

we obtain results shown in Fig. [3j Amazingly, the sign 
problem between the expansion orders is essentially elim- 
inated and accurate results can be obtained already for 
relatively small expansion orders ^ 5, making the entire 
scheme a viable solution to the regularization of diagram- 
matic expansions with zero-convergence radius (increas- 
ing N does not change the picture substantially). 

At this point we note that series convergence for 
fixed N allows us to truncate the series at some or- 
der p*(iV) such that the difference between /j 



The other important difference between Eqs. ( 15 1 and 



( 10 1 is that the y variables in Eq. ( 10 1 do not contain the 



J dx e ^Zp_^ j^(x) and /jv is smaller than \I^ 



P*,N = 
N-I{X)\. 

This would account for the regularization technique in 
the form Q when the set of approximations to the ex- 
ponential function is always in the form of a finite-order 
polynomial in x. Clearly, there is a lot of room for opti- 
mizing the form of such polynomials, and the success of 
the near-Gaussian resummation, which in essence per- 
forms such an optimization, proves the point — radical 
improvements in efficiency are possible along these lines 
even with relatively small p*. One of such approaches 
(tested for our model with quite satisfactory results) is 
based on polynomials that reproduce exactly the first 
moments of the distribution e"^""^^ . 

To deal with a generic non-local coupling between fields 
g-A|Vi|'|^2|' ^ g-XxiX2 propose to start with the fol- 
lowing regularization (later on it can be further modi- 
fied/optimized) 



-XxiX2 



lim 



9{yi)9{V2) 



\xiX2 

N 



f{Vl)f{V2 



(15) 

where yi — u'a:^^, and u' = (l/N"-)^/'^ with the exponent 
a < 1/2. The entire functions g{y) and f{y) have to sat- 
isfy criteria similar to those mentioned above: be even 
functions of y, decay on the real axis, and diverge not 
faster then e'^' in the complex plane. These requirements 
are satisfied if / and g are based on linear combinations of 
Bessel functions. We also demand that /(O) = g{0) = 1, 
the Taylor series expansion of g starts from an order y^'^ 
with as > 1, and the modulus of expression in paren- 
theses is smaller than unity. This also constitutes no 
problem using proper linear superpositions of J„. The 
resulting regularization is guaranteed to work and, most 
importantly, it can be generalized to deal with arbitrary 
interaction potentials (two-body, three-body, etc.) be- 
cause the / and g functions are defined in terms of local 
fields only. We had to proceed in this way instead of re- 
peating the solution of the single-field problem because 
the expansion in y^ cx |'0i||V'2| does not allow the dia- 
grammatic technique, contrary to yf cx IV'iP- 



parameter A, and the diagrammatic expansion goes thus 
in integer powers of A and u'^. We believe that with an 
appropriate choice of functions it is possible to proceed 
with a = 1/2 and u' = {X/Nf/^, just as before. This 
approach has the "advantage" of being based on a single 



parameter expansion. Likewise, one may apply ( 15 1 with 



a < 1/2 and an A- independent parameter u' to deal with 
the single-site regularization. 

In conclusion, diagrammatic series with zero conver- 
gence radius can be dealt with by approximating the 
original interaction exponential with an order-p polyno- 
mial in the integer powers of the original interaction and 
regularizing counter-terms. The action of the regular- 
ized theory is obtained by exponentiating the polynomial. 
The accuracy of the approximation is controlled by the 
parameter p. While the polynomial has to satisfy a num- 
ber of requirements guaranteeing, in particular, that the 
sequence of approximate theories approaches the origi- 
nal one at p — > oo (counterintuitively, this property re- 
quires special care) , there is still a continuum of possible 
choices, with some characteristic examples tested on the 
basis of zero-dimensional \'ip\'^ model. We believe that 
the outlined approach — especially in the context of Di- 
agrammatic Monte Carlo — opens up an opportunity to 
utilize Feynman's diagrams as a generic tool to address 
strongly correlated classical- and quantum-field systems. 
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